Pangenome-genotyped structural variation improves molecular phenotype mapping in cattle.

Expression and splicing quantitative trait loci (e/sQTL) are large contributors to phenotypic variability. Achieving sufficient statistical power for e/sQTL mapping requires large cohorts with both genotypes and molecular phenotypes, and so, the genomic variation is often called from short-read alignments, which are unable to comprehensively resolve structural variation. Here we build a pangenome from 16 HiFi haplotype-resolved cattle assemblies to identify small and structural variation and genotype them with PanGenie in 307 short-read samples. We find high (>90%) concordance of PanGenie-genotyped and DeepVariant-called small variation and confidently genotype close to 21 million small and 43,000 structural variants in the larger population. We validate 85% of these structural variants (with MAF > 0.1) directly with a subset of 25 short-read samples that also have medium coverage HiFi reads. We then conduct e/sQTL mapping with this comprehensive variant set in a subset of 117 cattle that have testis transcriptome data, and find 92 structural variants as causal candidates for eQTL and 73 for sQTL. We find that roughly half of the top associated structural variants affecting expression or splicing are transposable elements, such as SV-eQTL for STN1 and MYH7 and SV-sQTL for CEP89 and ASAH2 Extensive linkage disequilibrium between small and structural variation results in only 28 additional eQTL and 17 sQTL discovered when including SVs, although many top associated SVs are compelling candidates.

Adenylate kinase 9 is essential for sperm function and male fertility in mammals.

Despite passing routine laboratory tests for semen quality, bulls used in artificial insemination exhibit significant variation in fertility. Routine analysis of fertility data identified a dairy bull with extreme subfertility (10% pregnancy rate). To characterize the subfertility phenotype, a range of in vitro, in vivo, and molecular assays were carried out. Sperm from the subfertile bull exhibited reduced motility and severely reduced caffeine-induced hyperactivation compared to controls. Ability to penetrate the zona pellucida, cleavage rate, cleavage kinetics, and blastocyst yield after IVF or AI were significantly lower than in control bulls. Whole-genome sequencing from semen and RNA sequencing of testis tissue revealed a critical mutation in adenylate kinase 9 (AK9) that impaired splicing, leading to a premature termination codon and a severely truncated protein. Mice deficient in AK9 were generated to further investigate the function of the gene; knockout males were phenotypically indistinguishable from their wild-type littermates but produced immotile sperm that were incapable of normal fertilization. These sperm exhibited numerous abnormalities, including a low ATP concentration and reduced motility. RNA-seq analysis of their testis revealed differential gene expression of components of the axoneme and sperm flagellum as well as steroid metabolic processes. Sperm ultrastructural analysis showed a high percentage of sperm with abnormal flagella. Combined bovine and murine data indicate the essential metabolic role of AK9 in sperm motility and/or hyperactivation, which in turn affects sperm binding and penetration of the zona pellucida. Thus, AK9 has been found to be directly implicated in impaired male fertility in mammals.

Graph construction method impacts variation representation and analyses in a bovine super-pangenome.

BACKGROUND: Several models and algorithms have been proposed to build pangenomes from multiple input assemblies, but their impact on variant representation, and consequently downstream analyses, is largely unknown. RESULTS: We create multi-species super-pangenomes using pggb, cactus, and minigraph with the Bos taurus taurus reference sequence and eleven haplotype-resolved assemblies from taurine and indicine cattle, bison, yak, and gaur. We recover 221 k nonredundant structural variations (SVs) from the pangenomes, of which 135 k (61%) are common to all three. SVs derived from assembly-based calling show high agreement with the consensus calls from the pangenomes (96%), but validate only a small proportion of variations private to each graph. Pggb and cactus, which also incorporate base-level variation, have approximately 95% exact matches with assembly-derived small variant calls, which significantly improves the edit rate when realigning assemblies compared to minigraph. We use the three pangenomes to investigate 9566 variable number tandem repeats (VNTRs), finding 63% have identical predicted repeat counts in the three graphs, while minigraph can over or underestimate the count given its approximate coordinate system. We examine a highly variable VNTR locus and show that repeat unit copy number impacts the expression of proximal genes and non-coding RNA. CONCLUSIONS: Our findings indicate good consensus between the three pangenome methods but also show their individual strengths and weaknesses that need to be considered when analysing different types of variants from multiple input assemblies.

Complex histories of gene flow and a mitochondrial capture event in a nonsister pair of birds.

Hybridization, introgression, and reciprocal gene flow during speciation, specifically the generation of mitonuclear discordance, are increasingly observed as parts of the speciation process. Genomic approaches provide insight into where, when, and how adaptation operates during and after speciation and can measure historical and modern introgression. Whether adaptive or neutral in origin, hybridization can cause mitonuclear discordance by placing the mitochondrial genome of one species (or population) in the nuclear background of another species. The latter, introgressed species may eventually have its own mtDNA replaced or "captured" by other species across its entire geographical range. Intermediate stages in the capture process should be observable. Two nonsister species of Australasian monarch-flycatchers, Spectacled Monarch (Symposiachrus trivirgatus) mostly of Australia and Indonesia and Spot-winged Monarch (S. guttula) of New Guinea, present an opportunity to observe this process. We analysed thousands of single nucleotide polymorphisms (SNPs) derived from ultraconserved elements of all subspecies of both species. Mitochondrial DNA sequences of Australian populations of S. trivirgatus form two paraphyletic clades, one being sister to and presumably introgressed by S. guttula despite little nuclear signal of introgression. Population genetic analyses (e.g., tests for modern and historical gene flow and selection) support at least one historical gene flow event between S. guttula and Australian S. trivirgatus. We also uncovered introgression from the Maluku Islands subspecies of S. trivirgatus into an island population of S. guttula, resulting in apparent nuclear paraphyly. We find that neutral demographic processes, not adaptive introgression, are the most likely cause of these complex population histories. We suggest that a Pleistocene extinction of S. guttula from mainland Australia resulted from range expansion by S. trivirgatus.

Detecting turnover among complex communities using null models: a case study with sky-island haemosporidian parasites.

Turnover in species composition between sites, or beta diversity, is a critical component of species diversity that is typically influenced by geography, environment, and biotic interactions. Quantifying turnover is particularly challenging, however, in multi-host, multi-parasite assemblages where undersampling is unavoidable, resulting in inflated estimates of turnover and uncertainty about its spatial scale. We developed and implemented a framework using null models to test for community turnover in avian haemosporidian communities of three sky islands in the southwestern United States. We screened 776 birds for haemosporidian parasites from three genera (Parahaemoproteus, Plasmodium, and Leucocytozoon) by amplifying and sequencing a mitochondrial DNA barcode. We detected infections in 280 birds (36.1%), sequenced 357 infections, and found a total of 99 parasite haplotypes. When compared to communities simulated from a regional pool, we observed more unique, single-mountain haplotypes and fewer haplotypes shared among three mountain ranges than expected, indicating that haemosporidian communities differ to some degree among adjacent mountain ranges. These results were robust even after pruning datasets to include only identical sets of host species, and they were consistent for two of the three haemosporidian genera. The two more distant mountain ranges were more similar to each other than the one located centrally, suggesting that the differences we detected were due to stochastic colonization-extirpation dynamics. These results demonstrate that avian haemosporidian communities of temperate-zone forests differ on relatively fine spatial scales between adjacent sky islands. Null models are essential tools for testing the spatial scale of turnover in complex, undersampled, and poorly known systems.

Inter- and intra-archipelago dynamics of population structure and gene flow in a Polynesian bird.

Islands are separated by natural barriers that prevent gene flow between terrestrial populations and promote allopatric diversification. Birds in the South Pacific are an excellent model to explore the interplay between isolation and gene flow due to the region's numerous archipelagos and well-characterized avian communities. The wattled honeyeater complex (Foulehaio spp.) comprises three allopatric species that are widespread and common across Fiji, Tonga, Samoa, and Wallis and Futuna. Here, we explored patterns of diversification within and among these lineages using genomic and morphometric data. We found support for three clades of Foulehaio corresponding to three recognized species. Within F. carunculatus, population genetic analyses identified nine major lineages, most of which were composed of sub-lineages that aligned nearly perfectly to individual island populations. Despite genetic structure and great geographic distance between populations, we found low levels of gene flow between populations in adjacent archipelagos. Additionally, body size of F. carunculatus varied randomly with respect to evolutionary history (as Ernst Mayr predicted), but correlated negatively with island size, consistent with the island rule. Our findings support a hypothesis that widespread taxa can show population structure between immediately adjacent islands, and likely represent many independent lineages loosely connected by gene flow.

A test of island biogeographic theory applied to estimates of gene flow in a Fijian bird is largely consistent with neutral expectations.

Islands were key to the development of allopatric speciation theory because they are a natural laboratory of repeated barriers to gene flow caused by open water gaps. Despite their proclivity for promoting divergence, little empirical work has quantified the extent of gene flow among island populations. Following classic island biogeographic theory, two metrics of interest are relative island size and distance. Fiji presents an ideal system for studying these dynamics, with four main islands that form two large-small pairs. We sequenced thousands of ultraconserved elements (UCEs) of the Fiji bush-warbler Horornis ruficapilla, a passerine distributed on these four Fijian islands, and performed a demographic analysis to test hypotheses of the effects of island size and distance on rates of gene flow. Our demographic analysis inferred low levels of gene flow from each large island to its small counterpart and little or none in the opposite direction. The difference in the distance between these two island pairs manifested itself in lower levels of gene flow between more distant islands. Both findings are generally concordant with classic island biogeography. The amount of reduction in gene flow based on distance was consistent with predictions from island biogeographic equations, while the reduction from small to large islands was possibly greater than expected. These findings offer a hypothesis and framework to guide future study of interisland gene flow in archipelagos as the study of island biogeography progresses into the genomic era.